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PREFACE 


The  Hardness  Verification/Surveillance  Testing  and  Standards  Program  (HAVSTAP) 
is  funded  by  the  Defense  Nuclear  Agency  (DNA)  and  other  government  agencies 
through  the  Corps  of  Engineers  (CoE)  U.S.  Army  Engineer  Division  Huntsville  (US- 
AEDH).  The  contract  monitor  for  the  CoE  is  Mr.  John  Loyd,  and  the  DNA  technica 
monitor  is  Mr.  Sam  Mauch. 

Mission  Research  Corporation  (MRC)  is  the  prime  contractor  for  the  HAVSTAP, 
and  IRT  Corporation  and  BoozAllen  and  Hamilton  Inc.  ar^^ontractors  to  MRC. 
Personnel  from  all  three  companies  were  instrumental  in  tl^^^^^^test. 
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The  Harry  Diamond  Laboratories  (HDL)  is  conducting  a  HEMP  hardness 
verification  prograi^Bi^As  part  of  thi.  program,  MH^STD-285  shielding 

effectiveness  tests  were  performed,  and  the  conducted  penetrations  _ _ 

H^^ere  direct  driven  with  high  level  pulses.  For  more  information 

on  the  HDL  effort  see  Reference  2. 

To  augment  the  HDL  tests ,  the  Defense  Nuclear  Agency  (DNA)  sponsored 
a  low  level  CW  test^^g^g^Bi  The  CW  test  was  conducted  from  26  September 
to  21  October  1988  to  determine  the  HEMP  induced  stresses.  The  measurements 
were  acquired  with  the  DNA  CWMS  HI  (Continuous  Wave  Measurement  System  HI). 
This  test  was  performed  under  the  auspices  of  the  Hardness  Verification/Surveillance 


Testing  Mid  Standards  Program  (HAVSTAP).  The  test  execution  mtd  lessons  learned 
Le  documented  in  the  Test  Director’s  Report  (Reference  1).  The  purpose  of 
report  is  to  describe  the  data  processing  and  to  analyze  the  results. 


1.2  TEST  OBJECTIVES 


The  specific  objectives  of  the  CWMS  tes 

1.  Acquire  HEMP  stress  estimates  at  internal  equipment  interfaces. 


2.  Determine  the  dominant  source  for  these  stresses. 

3.  Complement  the  DNA  HEMP  stress  date  base  for  Ground  Based  C8I  facilities 
with  stress  data  from  a  well  shielded  facility. 


4.  Demonstrate  the  utility  of  CWMS  testing  for  well  shielded  facilities. 


6.  Provide  recommendations  for  additional  testing  to  support  the  HDL  hardness 
verification  statement. 


To  meet  these  objectives,  C W  illumination  and  CW  direct  drive  tests  on  the 
conducted  penetrations  were  performed.  Table  1  provides  an  overview  of  the  number 
of  test  points  monitored  during  the  various  elements  of  the  test.  CW  illumination  is 
discussed  in  Section  2,  the  single  line  direct  drive  data  are  described  in  Section  3,  and 
bulk  direct  drive  results  are  contained  in  Reference  3. 
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Table  !.  Number  of  teet  location,  for  Illumination,  .ingle  line  direct  drive, 
and  bulk  direct  drive  tests. 


(b)  Valuable  lessons  were  learned  in  HEMP  testing  and  analysis. 
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SECTION  2 

ILLUMINATION  MEASUREMENTS 


2.1  EXTRAPOLATION 


The  illumination  test  was  performed  with  the  CWMS  inverted-V  antenna. 
Due  to  the  uneven  terrain  and  the  security  fence  around  the  property  there  were  only 
few  options  for  installing  the  antenna.  Figure  2  shows  the  antenna  location  for  the 
CW  illumination  test  at  the  ETC.  The  inverted-V  antenna  simulates  a  horizontally 
polarized  incident  plane  wave  with  elevation  angle  roughly  34°  with  field  components 

E„,  H„  and  H,. 


Figure  2  also  shows  the  location  of  the  reference  sensor  (on  the  antenna 
center  line,  20  m  south  of  the  antenna,  1  m  above  ground).  The  reference  sensor 
measures  the  total  H,  field  produced  by  the  radiating  antenna.  The  reason  for  this 
particular  placement  of  the  reference  sensor  is  to  minimize  undesirable  scattering  off 
the  building.  Thus  the  reference  sensor  is  placed  on  the  opposite  side  of  the  antenna 
utilizing  the  field  symmetry  (with  respect  to  the  y-axis)  properties  of  the  antenna. 


The  field  measured  by  the  referenc^en^^s  identical  to  the  field  which 
would  be  at  the  symmetrical  location  (i.e.,  t^ie  absence  of  the  building. 

CWMS  HI  measures  the  current  induced  on  cables,  relative  to  the  simulation 
field  at  the  reference  sensor.  Mathematically  each  measurement  is  a  transfer  function 
of  the  form: 


A 

Ajm 


(1) 


where  TP  represents  the  test  signal  induced  on  the  cable  measured,  HZ(R)  is  the  x- 
component  of  the  magnetic  field  measured  at  the  reference  sensor,  and  the  subscript 
“ill”  denotes  an  illumination  measurement. 

To  obtain  an  estimate  of  the  HEMP  induced  current,  the  transfer  function 
must  be  multiplied  by  the  total  (incident  plus  ground  reflected)  HEMP  field  at  the 
reference  sensor  location.  Specifically,  this  field  is  the  total  Hs  which  would  exist  at  the 


r34 ^he  Sn“!f  pLSnlt^elative  permittivity 
e/eo  =  10,  and  conductivity  tr  =  4  nunho/m. 

Since  the  simulation  fields  satisfy  Maxwell's  equations,  the  other  field  com¬ 
ponents  are  Somatically  scaled  to  the  appropriate  HEMP  euvronment. 

The  incident  HEMP  waveform  is  the  double  exponential  waveform 


(V/m) 


(2) 


£y~(t)  =  s>(e“"  -  «■*) 

V  .k,  (t,.  values  for  the  parameters  are  provided.  The  relationship 
l"  IhHavrfo^Twith  other  pMP  environments  is  discussed  in  Reference  4.  The 
corresponding  magnetic  field  is 


jyqne  _  pne JtfQJ  (A/m) 

and  the  magnetic  field  component  in  the  x-direction  is 

ginc  _  ffinceos( 34°)  U/m). 

EARLY  TIME  HEMP 


(3) 


(4) 
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Figure  4  shew.  H~.  the  correepontog  Wtel  H“.  which  wee  used  to 
extrapolate  all  illumination  test  measurements  to  HEMP. 

Butterworth  Band  Pass  Filter:  The  extrapolation  is  performed  fat  the  fr* 
ouency  domain.  To  obtain  time  domain  transients,  the  frequency  spectra  must  be 
favei2  Fourier  transformed.  CWMS  illumination  measurements  are  acquired  from 
100  kHz  to  100  MHz,  thus  the  extrapolated  spectra  are  band  limited  to  this  range. 
This  limitation  is  imposed  by  the  usable  frequency  range  of  the  simulator  antenna, 
and  a  justification  for  this  frequency  range  is  beyond  the  *cope  of  this  report  (e.g., 
Reference  5).  The  practical  consequence  of  this  limitation  is  that  ^ 

Fourier  transforms  can  suffer  due  to  the  abrupt  transition  of  the  spectra  to  zero  £100 
kHz  and  100  MHz.  To  smooth  the  transition  at  the  ends  of  the  measurement  band, 
a  Butterworth  FUter  (pass  band  200  kHz  to  80  MHz)  is  used  on  all  cable  current 

measurements.  (Therefore,  Fourier  transforms  5  ^ 

of  the  Butterworth  Filter  characteristic  with  the  HEMP  H£*  from  Figure  4  is  s 

in  Figure  5. 

Thus  the  formula  for  extrapolating  illumination  measurements  is: 


/<*'(/)  '  <5) 

where  TP  indicates  the  test  point  measured,  H,(R)  is  the  simulation  field  measured 
by  the  reference  sensor,  and  is  as  in  Figure  5. 

Simulation  Error.  The  incident  HEMP  is  specified  as  a  plane  wave.  How¬ 
ever,  the  inverted-V  antenna  can  only  approximate  a  plane  wave.  The  extrapolation 
numerically  simulates  exactly  the  specified  HEMP  environment  at  the  reference  sensor, 
and  also  at  the  corresponding  locatio  JgBH^^due  to  the  field  symmetry  with 
respect  to  the  y-axis).  At  all  other  locations,  the  simulated  fields  only  approximate  the 
specified  HEMP  environment.  The  accuracy  of  the  approximation  was  mvestigated  m 
References  6  and  7.  _ .  _ — — 


AMPS/M/H* 


Figure  4.  Incident  HEMP  Hy*  vs.  total  HEMP  H‘/‘  spectra,  and  total 
HEMP  H‘0<  time  waveform. 
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AMPS/M/Hz 


Figure  5.  Extrapolation  function  including  Butterworth  Filter  for 

cable  current  transfer  functions. 

2.2  EXTERNAL  MEASUREMENTS 

2.2.1  EXTERNAL  FIELD  MAP 

The  purpose  of  the  field  map  at  the  external  locations  FM1,  FM^an^j^ 
is  to  provide  data  on  the  non-uniformity  of  the  simulation  field  along  th^^^^^F** 
MQP'CWMS  measures  field  transfer  functions  analogous  to  cable  current  transfer 
Unctions.  In  other  words,  just  like  cable  measurements,  fields  are  also  measured 
relative  to  the  H,(R)  measured  at  the  reference  sensor.  Therefore,  to  extrapolate  field 
transfer  functions  to  HEMP,  the  measurements  are  multiplied  by  HEMP  Hx°  shown 
in  Figure  4.  (The  Butterworth  Filter  is  not  necessary  for  field  measurements  because 
their  low  frequency  content  is  small,  and  at  the  high  frequency  end  the  rolloff  of  the 
HEMP  spectrum  automatically  provides  a  sufficiently  smooth  transition.) 

Figure  6  shows  the  extrapolated  field  measurements.  Ideally  the  field  at 
the  three  field  map  locations  would  be  identical  to  the  field  transient  at  the  reference 
sensor  shown  in  Figure  4. 
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to  time  (t  =  0  when  the  field  arrives  at  the  reference  sensor)  due  to  the  differences  m 
propagation  length  to  the  three  locations.  (The  high  frequency  components  emanate 
from  the  bicone  at  the  center  of  the  inverted-V  antenna.) 


Consequently  the  currents  induced  on  thJIfrin  the  illumination  test  are 
not  equal  to  those  which  would  be  induced  by  a  HEMF  plane  wave.  To  bound  this 
illumination  deficiency,  the  field  measurementsg^m^re  averaged  as  follows: 


AVE  is  the  average  simulation  H,  along  th^PJH^p  relative  to  the  H, 
at  the  reference  sensor,  i.e.,  AVE  is  the  average  illumination  deficiency.  AVE  would 
be  equal  to  unity  at  all  frequencies  if  the  simulated  field  were  a  uniform  plane  wave. 
The  reciprocal  of  AVE  quantifies  the  simulation  field  deficiencies,  i.e.,  the  correction 
function  to  compensate  for  the  deficiencies  is  given  by: 


CORR(f)  =  1/AVE(/) 


(7) 


The  magnitude  of  this  correction  function  is  shown  m  Figure  7.  Note  that 
at  frequencies  below  1  MHz  the  deficiency  is  only  on  the  order  of  4  dB,  but  above  10 
uL  the  error  is  greater  than  10  dB  and  rapidly  varying  with  frequency.  The  reason 
for  the  spiky  behavior  is  that  AVE  is  the  sum  of  three  complex  spertra  which  are 
Fourier  transforms  of  delayed  transients  (cf.  Figure  6).  The  effect  of  the  tone  delays 
between  the  three  transients  is  that  at  every  frequency  their  phases  either  combine 
constructively  (creating  a  peak),  or  destructively  (producing  a  null). 

XI  TEST 

Hx  (REF) /Hx (AVE)  _ . 


2.2.2 

This  section  discu 
illumination  test. 


the  currents  measured  on  th 
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The  cable  cnnent  meagorementa  weremadeon  thedlrty  eid.of  U>«buUdmt 
AWd  (in  the  doghou.es).  The  waveform  nonn.  for  th«  current.  ere  contented  » 


A _ JU,  A 


75 


u  rtrt 

ILLUKHUTIOW  PATA 


17 


18 


MAGNITUDE  dB  (A/H») 
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MAGNITUDE  dB  IA/H*) 
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Figure  12  shows  the  total  current  The  currents  induced 

~  iighe4Hlf  than  expected.  Although  the  currents  measured  on  the 
are  probably  not  very  accurate  estimates  of  the  actual  currents  induced  by 
a  HEMP,  there  are  no  data  to  bound  the  error.  This  is  further  discussed  in  Section 

2.3.I.3. 


2.3  INTERNAL  MEASUREMENTS 


The  illumination  measurements  _di 
rents  inside 


RTilVTliiF 


In  addition,  internal  field  measurements  were  taken  at  two  location 


of  the  illumination  cable  current  measurements  is  to  dete 


These  inte 


he: 


MAGNITUDE  dB  (A/Hz) 


K1  TEST 
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TestcmnlsJabeled  “B”  are 

and  test  points  labeled^”^^|^HH^^^I^^P 
“^H^^KKiielded  and  unshielded  cables  and  wire  bundles.  Current  mea- 
P°m  *  -1  wires  are  generally  impractical  at  operating  facilities  because  they 

SKSSKi-  L  CWMS  test  signals  are  usually  dominated  by 
Uie  operating  signal,  which  constitutes  noise  to  CWMS  measurements.  Moreover  .t 
is  impossible  to  survey  complex  communications  facilities  on  t  e  asis  o  sing  e 
currents. 


In  the  following  the  quality  of  the  measurem 
stress  estimates  are  presented,  and  the 
Section  2.2.2  is  estimated 


discussed,  the  HEMP 


illumination  deficienc 


2 .3.1.1  D ATA  QXJ ALITY  S/N Ratio:  Since  the  ETC  was  operating  throughout 
the  CWMS  test,  the  operating  currents  on  the  cables  measured  were  competing  with 
the  CWMS  induced  test  signal.  Therefore,  at  every  internal  test  point  the  facility 
generated  ambient  noise  was  measured  (with’ the  CWMS  antenna  nmtad),  and  com- 
pared  with  the  CWMS  signal  at  the  test  point.  This  Signal-to-Noise  {S/N)  relation  is 
a  comparison  of  the  current  detected  in  the  signal  (S)  channel  of  the  network  analyzer 
due  to  CWMS  illumination,  and  then  due  to  the  ambient  noise  (N).  Note  that  S  is  no 
a  CWMS  transfer  function,  but  is  the  numerator  of  the  transfer  function  (the  units  of 

S  and  N  are  A/y/Hz). 

ical  examples  of  the  signals  and  noise  at  test  point  ^ _ 

'  (Strictly  speaking,  the  quantities  compared  are  noise, 
ana  signal  plus  noise  because  the  noise  is  always  present.)  The  Signal-  to-Noise  ratio 
(in  dB)  is  obtained  by  subtracting  the  dotted  line  (Noise  in  dB)  from  the  solid  line 
(Signal  in  dB)  in  Figure  13.  However,  instead  of  plotting  this  S/N  ratio  vs.  frequency, 
it  more  illustrative  to  compare  S  and  N  directly  as  shown.  With  this  illustration,  for 
example,  it  can  be  observed  how  the  noise  “bites”  the  signal  at  isolated  frequencies, 
and  thus  creates  spikes  in  the  transfer  function. 
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dB(A/\//7»)  dB(A/y/H») 


2 


-  -  *  **  «  " 

Y/l Th^forc,  the  following  figure  of  merit  wan  developed  for 
S/N  quantification: 


riVSN  -  5J  £  /(") 

nssl 


(dB) 


where 


/(n)  =  S/N  if  S/N  >  6dB 
-  0  if  S/N  <  6dB 


[«irrmvg 


^Fourier  Transforms:  If  the  signal  to  noise  ratio  is  low  for  a  mea- 
sorement,  the  true  (but  not  measurable)  signal  lie.  below  the  moused  nowe  Thus, 
the  measurement  constitute.  »n  upper  bound  on  the  true  signal  at  the  test  point. 
Noisy  measurements  can  .till  be  extrapolated  and  Fourier  transformed,  but  the  re- 
suiting  transient  may  not  be  accurate.  A  measure  of  this  inaccuracy  is  obtained  by 
performing  sine  (or  odd)  and  cosine  (or  even)  transforms,  and  comparing  the  resulting 
transient.  Ideally  they  would  be  equal.  This  is  illustrated  in  Figure  15,  which  shows 
the  even  and  odd  transforms  for  two  extrapolated  currents  measurements. 


Again  a  figure  of  merit  is  needed  to  perform  an  automated  sift  through  large 
data  bases  for  Fourier  tranform  quality  control.  Such  a  figure  of  merit  is  a  Fourier 
tranforms  error  (FTE)  defined  by: 


1  ft  (tm)  ~  /o(*m) 


i 


PKI 
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KO.  - - 

average  signal  to  noise  ratio 

15  I  ■  ■  '  1  ill'll  ' 


12 


FTE  is  a  unitless  quantity  and  measures  the  area  between  the  even  and 
odd  transforms,  normalized  to  the  peak  of  the  odd  transform.  Therefore,  perfect 
transforms  have  FTE  =  0.  As  a  rule  of  thumb  based  on  judgment,  good  transforms 
have  FTE  values  less  than  0.3* 


One  would  expect  a  causal  relationship  between  Signal-to-Noise  ratios  and 
quality  of  Fourier  transform.  Specifically,  for  poor  S/N  (small  AVSN)  the  quality  of 
the  transform  should  also  be  poor  (large  FTE).  The  scatter  plot  of  FTE  vs.  AVSN  m 
Figure  16  confirms  that  the  two  figures  of  merit  are  correlated. 


Upper  Bound:  Upper  bounds  exist  for  several  waveform  norms,  which  can 
be  used  to  bound  poor  Fourier  transforms.  For  example,  the  peak  is  bounded  by  the 
norm  PINT: 


PKI  =  max  |  fo{t)  |=  -  |  ffl  F{<jj)Jwidu  |<  —  f  |  -F(u>)  |  du>  -  PINT  (A)  (11) 

t  7T  J fo  Jfo 
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OxlO"3  * 

OxlO-* 


200x10-* 


400x10 


600x10-* 


600x10- 


1000x10- 


TIME  (SEC) 
EXTRAPOLATED  TO  HEMP 


XI  TEST 


TP:  B66 


Figure  15.  Even  and  odd  tranform  examples  for  a  good  transform  (TP: 

B63,  FTE  =  0.05),  and  the  poorest  tranform  (TP:  B86,  FTE 
=  1.54).  (FTE  calculated  per  equation  (10)). 


29 


►was 


not  illuminated  with  a  plane  wave. 


3.  T 


mlation  anteima  was  only  located  in  one  position,  namely  to  the  south 
^^^^JUPoriented  east-west.  Therefore  the  building  and  the  penetrations  on 
thenorthside  were  on  the  shadow  side,  i.e.,  not  strongly  illuminated. 

4.  The  measurement  quality  is  not  perfect  at  all  internal  test  points  due  to  the 
presence  of  ambient  noise. 

These  error  sources  are  discussed  in  the  following. 

was  shown  in  Sections  2.2.2.1  and  2.2.2.2  that  the  extrap- 
olated  extemaj^pcurrents  are  13  to  16  dB  too  low.  Consequently,  the  HEM^ 
current  estimatesmay  be  low  at  internal  test  points  which  are  connected 

Since  the  correction  factor  CORR(f)  defined  in  Section  2.2.1 
BSSflSstSnHKh  frequency  but  varies  considerably,  a  more  accurate  estimate 
of  the  error  is  obtained  by  applying  CORR(f)  to  internal  measurements. 

Table  2  lists  the  peak  HEMP  current  estimates  (from  Section  2.3.1.2)£ 

The  spectra  of  these  11  currents  were  then  multiplied  by 
product  inverse  transformed.  The  column  entitled  “Scaled 
by  CORR(f)"  contains  the  resulting  peaks.  The  last  column  shows  the  ratios 
expressed  in  dB.  Thus  the  average  error  is  11.7  dB  or  a  factor  of  4.  Therefore  the 
error  of  the  HEMP  stress  estimates  at  East  SAMT  test  points  is  approximately 
12  dB.  This  bound  is  conservative  because  CORR(f)  overestimates  the  effect  of 
the  illumination  deficiency  (cf.  Section  2.2.1).  This  is  the  only  appHc^^rror 
because  the  three  SAMTs  are  electrically  well  isolated  as  shown  ^^Vbulk 

Since  the  inverted-V  antenna  used  in  the  illumination  test  has  been  shown  to  be 
symmetric  with  respectt^h^eM^^^^ds),  the  12  dB  error  bound  also 
approximately 

^situation  becomes  more  complicated  because 
these  test  points  are  connected path,  not 
hardwire)  ft 

certainly  be  conservativ^^^ppl^^^^d^erro^Dound  to  all  stress  estimates  in 
these  two  rooms.  However,  it  would  nrahahhi  also  ^«>  mitt»  trnrpahstic  because 
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2  £||^BHFwas  dearly  not  illuminated  with  a  plane  wave,  but  the  currents 
induced  are  not  sm 


be  necessary  to  compare  the  UW  mummauon  measur«u«i»  - - — - — 

and  reliable  data.  However,  for  the  particular  HEMP  illumination  simulated 

it  is  expected  that  this  error 

SS^K^SaU^n^erhap^veiwiegative,  indicating  that  the  CW 
illumination  results  are  actually  overestimates. 


3.  The  CW  illumination  did  not  adequately  test  the  integrity  of  the  north  building 
wall  because  only  a  single  antenna  position  was  possible  within  the  time  con¬ 
straints  of  the  test.  (Data  on  the  shielding  effectiveness  of  the  north  wall  were 
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the 


^^^^^^^—^hereforeTthe  stress  contributions  through 
^^S^mauTbut  in  any  event  cannot  be  quantified  with  the 
alone  except  for  the  telephone  penetration. 


test  results 


4.  Ambient  Noise:  As  discussed  in  Section  2.3.1.1,  the  facility  generated  no.se 
interfered  with  the  CW  test  signal  However,  since  the  measurements  are  signa 
plus  noise,  the  stress  estimates  are  automatically  conservative,  and  no  additional 

error  bound  is  needed. 


For  measurements  yielding  poor  Fourier  transforms,  the  norm  PINT  is  used 
to  bound  the  peak  current.  Again,  PINT  is  already  an  upper  bound,  and  thus  no 
additional  error  bound  is  required. 


In  summary,  12  dB  (factor  of  4)  is  a  consOTgta^mjra^gjjj^for  t^ie  HEMP 
stress  estimates  in  Section  2.3.1.2,  except  perhaps  ^^^^^^■■Fhere  the  error 
cannot  be  defined.  Thus,  the  histogram  in  Figure™  would  be  shifted  by  a  factor  of 
four  to  bound  the  errors. 


2.3.2  INTERNAL  FIELD  MEASUREMENTS 


During  the  CW 
tions  (Dl,  D2 


internal  fields 


measured  at  two  loca-l 
I  were  taken  1  m 


above  the  floor 


A  measure  of  the  global  shielding  effectiveness  is  obtained  by  comparing 
the  internal  fields  with  the  external  fields.  The  dominant  components  of  the  incident 
field  are  E„,  Hx,  and  H,.  The  fields  inside  the  building  are  produced  by  diffraction, 
scattering,  and  reradiation  from  cable  currents,  hence  there  are  no  a  priori  preferred 
or  theoretically  dominant  field  components.  Therefore,  for  completeness  one  would 
have  to  measure  all  six  field  components.  However,  such  an  exhaustive  investigation 
into  the  merits  of  internal  field  mapping  was  not  possible  within  the  test  schedule. 


f *1 

i+v 


^  3i  AfUr, 
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There  are  important  open  technology  questions  associated  with  internal  field 
mapping,  such  as: 

.  Is  the  CWMS  capable  of  verifying  compliance  with  a  100  dB  shielding  require¬ 
ment?  (SpecificX,  »  range  ?d  noise  raject.™  features  suffice* 

to  measure  fields  100  dB  below  the  external  simulation  fields?) 

•  How  many  internal  locations  must  be  mapped,  and  where? 

•  How  well  are  internal  field  map  measurements  correlated  with  MIL-STD-  285 
measurements? 

The  internal  field  map  ™  Perform'<i  °"'y  °"  “ 

experimental  basis  to  determine  the  usefulness  of  such  measurements.  Therefore, 
OT*y  the  components  identified  in  Table  3  were  measured.  Note  that  there  are  three 
configurations  tor^jH/^ooi  AS  (closed,  open,  and  taped).  Only  measuremaits 
acquired  in  the  noSKonfiguration  (door  closed)  are  drscussed  m  this  section,  the 
other  two  axe  deferred  to  the  next  section. 


Data  Quality:  The  ambient  noise  discussed  in  Section  2.3.1.1  also  affects 

*  J  _ .  .a  %  aa  _ ire 


(PINT) 

a Wjj^M J^SmpTete  list  is  contained  in  Table  4,  Section  2.3.3. 


Table  3. 


totemal  field  component,  meamred  during  CW  illumination 
testing. 


"Internal  Field  Field  Components 

Man  Location  Configuration  Measured - 


D1 

D1 


D1 


n/a 


E**  E«  Eg,  H* 
E„E. 

E*,  E* 

E„  E„  H,,  H* 
E, 

E«,  E, 

Ex,  E,,  H«,  H, 


Field  Attenuation  (FATT):  A  measure  of  shielding  effectiveness  or  field  at, 
tenuation  is  obtained  by  dividing  the  internal  field  measurements  by  the  ®  • 

It  seems  that  the  most  unambiguous  form  of  the  external  field  is  e  inci 
Thus  the  attenuation  of  electric  fields  (EFATT)  b  defined  by: 


EaFATT{f)  =  -20  log 


Eine 


(dB) 


(12) 


where  E’."*  ‘is  one  of  the.  three  internal  electric  field  components.  Similarly,  the  atten 
uation  of  magnetic  fields  (HFATT)  b  defined  by: 


H.FATTU)  =  -20 


[dB) 


(13) 


where  Eint  b  one  of  the  three  internal  magnetic  field  components.  Note  that  these  def¬ 
initions  are  independent  of  ground  conductivity,  angle  of  incidence  (antenna  standoff 
dbtance),  and  other  test  conditions. 
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FATT  is  to  divide  the  extrapolated  measured  by  the  incident  HEMP  spectrum. 
For  example,  EFATT  for  the  x-component  is 


E.FATT(f)  =  —20  log 


/  > 

n? 

[b.W, 

(dS) 


B,FATT(f )  =  -20  log 


if  B?  \ 


376.7 

£ine 


(dB) 


The  FATT  was  calculated  for  all  four  locationqflH0^|^P|HpHHREr 
the  results  are  shown  in  Figures  21  to  24.  *or  comparison,  me  FATT 
TJJSSSs  of  MIL-STD-188-125  (Draft)  axe  included. 

These  figures  show  that  the  100  dB  electric  field  attenuation  cannot  be 
demonstrated^HI^MJbecause  the  internal  field  measurementsu^iois^om- 
inated  as  seen  from  Figures  19  and  20^  Moreover,  at  fre 
’da 


ment  was  met 

ic  field  attenuation  required  by  MIL-STD-188- 


Thus,  the  CWMS  data  acquisition  equipment  possesses  a  sufficiently  great 

^ynamif  range  to  demonstrate  compIiance_with_the_pro£osed_fieId_attenuation  curves 

In  a 

^visioiM>ntfIIr^TD-188-12^h^^^^^Si  attenuation  requirement  was  limited  to 
the  range  5  MHz  to  100  MHz.  fl|HBPhe  CWMS  was  able  to  demonstrate  100 
dB  electric  field  attenuation  ove^tusmodified  range  despite  the  noise. 
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ATTENUATION  dB  ATTENUATION  dB 


FREQUENCY  (H*) 

SOLID:  Hs  FATT  DOTTED:  H,  FATT 
DASHED:  MIL-STD-1M-1J5  (DRAFT) 


Figure  21.  Electric  (E.FATT,  F  FATT j  an^naglietic  field  attenuation 
(HrFATT,  H,FATT)^B^^HBRlocation  Dl, 

vs.  FATT  required  per  MIL-STD-188-125  (Draft). 
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ATTENUATION  dB 


FREQUENCY  (H*) 

SOLID:  H,  FATT  DOTTED:  H,  FATT 
DASHED:  MIL-STD-IM-KS  (DRAFT) 


Figure  22. 


Electric  (E,FATT,  E^^TTWn^gametic  field  attei^^toi 

Wrr,  HsFATT)^BHIHBV  (location  D2,4^MI 
[  vs.  FATT  required  per  MIL-STI)-188-125  (Draft). 
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ap  NoiivnNaiiv 


ATTENUATION  dB 


FREQUENCY  (Hi) 

SOLID:  E,  FACT*  DDTTED:  E,  FATT 
niCUPTi-  MIL-STD-lM-125  (DRAFT) 


FIELD  ATTENUATION 


FREQUENCY  (H*) 

SOLID:  H,  FATT  DOTTED:  H.  FATT 
DASHED:  MIL-STD-1M-1JS  (DRAFT) 


Figure  24. 


Electric  (E.FATT,  E,FATT)  and  magnetic  field  attenuation 
(H-FATT,  H,FATT)4HHH^H^  (location  Wl,^|| 
j  yg.  FATT  required  per  MXL-STD-188-125  (Draft). 


faej 

in 


c  A&'l 

1 1*3  i  •"e-'hj . 
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K1  TEST 


FREQUENCY  (Hr) 
TP:  02 

UNPROCESSED  DATA 


Figure  26. 


Internal  fields  Ey**  at  test  point  D2; 
(unextrapolated  measurements). 
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ATTENUATION  dB  ATTENUATION  dB 


FIELD  ATTEHUAiiOn 


FIELD  ATTENUATION 


Figure  27.  Elec^ic  field  attenuation  at  locations  D1  and  D2 

mm  vs*  required  per  MIL-STD-188-125  (Draft). 


f'ffc jc.  H6>  de(c.f<r*{ 

i  a  i't'S  e  t'Cf~ i  *j 
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FREQUENCY  (Hz) 
TP:  C31 

UNPROCESSED  DATA 


XI  TEST 


FREQUENCY  (Hz) 
TP:  C32 

UNPROCESSED  DATA 


Figure  28.  Unextrapolated  measurements  of  cable  currents  at  test  points 
G31  and  G32  for  various^^mH^bonfigurations.  The 
measurements  are  plotted  from  10  MHz  to  100  MHz  to 
enhance  legibility.  The  narrow  spikes  are  due  to  ambient 


noise. 
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NUMBER  NUMBER 


70 
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40  ■; 
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20  •; 
10 

0 


■■■I1- 
1.5  2.0 

PINT  (AMPS) 
SURVEILLANCE  TEST  2 
ILLUMINATION  DATA  (INTERNAL) 


2.5 


3.0 


N  ■  62 

MEAN  s  34.11  bA 

—l — I 1 r— 
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PINT  (AMPS) 

K1  TEST 

ILLUMINATION  DATA  (INTERNAL) 


Figure  29. 


Internal  HEMP  stresses 

terms  of  PINT  (both  populations  extrapolated 

to  HEMP  defined  in  equation  (15). 
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NUMBER 


mean  *  1.66E-4 
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Since  the  bulk  direct  drive  technique  is  still 
■■^■SltagerS^rsingle  line  direct  drive  results  are  discussed  here; 
the  bulkdirect  drive  data  are  evaluated  in  a  subsequent  report  (Reference  3). 


Single  line  direct  drive  measurements  are  complex  transfer  functions  of  the 


form: 


™=(£) 


SDD 


(A/A) 


(16) 


where  the  numerator  (signal)  is  the  current  measured  at  the  internal  test  point  (TT), 
and  the  denominator  is  the  current  driven  on  the  shield  of  the  extwnal  cable  Ai. 
subscript  “SDD"  indicates  that  the  measurement  is  a  single  line  direct  drive.  Sine 
lie  transfer  function  was  obtained  using  low  level  CW,  equation  (16)  only  predicts  the 
linear  transfer,  i.e.,  it  does  not  account  for  nonlinear  protection  features. 


3.1.1  INSERTION  LOSS 

The  measured  transfer  functions  ar^h^ttenu^io^rovide^j^h^ene- 

tration  hardening  as  a  function  of  frequency.^HHHHU^Rl^HHpiH^HHmmHIB’ 

►Figure  31  shows  two  typical  examples. 


Note  that  the  attenuation  varies  considerably  with  frequency.  Therefore,  the 
attenuation  cannot  be  readily  andcompletely  described  with  a  single  numerical  value. 
To  comprehend  and  compare  the®measurements  it  is  necessary  to  condense  the  data, 
i.e.,  another  figure  of  merit  is  needed. 

Since  there  does  not  appear  to  be  a  standard  method  of  achieving  this,  after 
considering  the  options,  the  following  figure  of  merit  for  insertion  loss  is  defined: 

IL[TP/Ai)  =  —20  l°9j~fQJjl  I  I  df  (17) 


which  indicates  the  ex 
where  fo  and  fi 


^1  test  point  (TP);  and 
Mathematically  IL 
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b  Simply  a  lineal  average.  Physically  IL  is  the  amount  by  which  the  penetration 
protection  would  attenuate  a  current  spike  (delta  function)  on  the  external  cable. 


icate  that  the  internal  test  point  is  hardwire 
example,  test  points  Bll  and  A4  are  the  sam 
and  A4  is  outside. 


In  addition  to  these  hardwi 


Values  printed  in 


ut  Bll  is  insid 


t  points,  a  number_of_cross 


Most,  but  not  al^jjtt^eld  currents  are  shunted  to  ground  at  the  penetration 
panel.  Therefore  part  of  ^^^Hmrrent^re  conducted  to  internal  equipment.  This 
part  is  the  land  line  dri 

To  determine  the  land  line  drive,  tte  direct  drive  measurements  must  be 
scaled  to  the  actual  currents  induced  on  th^^cable^^hyHu|^^^^^^  Thus 
the  land  line  drive  (LLD(TP))  at  an  internal  test  point^^^PH^^HBV  8iven 

by: 

{A/Hz)  (18) 

The  various  entries  in  expression  (18)  are  defined  by  equations  (16)  and  (5)  in  Section 

2.1. 

Physically  equation  (18)  adds  the  land  line  drive  contributions  from  each  of 
cables^These  contributions  are  scaled  to  precisely  the  actual  currents 
measured  on  th^m^t  cables  in  the  illumination  test  (cf.  Section  2.2.2. 1). 

Therefore,  the  land  line  drive  should  theoretically  always  be  smaller  than  the 
HEMP  stress  obtained  from  the  illumination  test  because  it  includes  land  line  drive. 
(Exceptions  could  occur  if  the  direct  drive  test  somehow  excited  different  coupling  paths 
not  driven  in  the  illumination,  but  this  does  not  apply  here.)  The  difference  between 
the  HEMP  stress  and  the  land  line  drive  is  the  local  coupling. 


S')  Ae'le.'fed. 


i  r\  i  e  f  i| 
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90  i 


„  the  land  line  drive  calculated  per  equation  U8)L  — 
compared  with  the  HEMP  stresses.  Table  7  contains 
■"KSSsforthe  HEMP  stresses  and  land  line  drives  shown  in  Figures 


32  to  34. 


It  is  evident  that  the  calculated  land  line  drives  exceed  the  HEMP  stresses, 
contrary  to  expectations.  The  reason  appears  to  be  that  the  measurements  are  imper¬ 
fect:  Note  that* calculation  of  the  land  line  drive  at  each  test  point  involves  9  direct 
drive  and  9  illumination  measurements,  none  of  them  perfect.  It  is  no  unreasona 
to  expect  that  the  cumulative  error  would  amount  to  the  order  of  10  <lg- 


3.1.3 


l<TER  INSERTION  LOSS 


Since  was  not  operational,  the  power  filte^va^^gs^^^^ 

direct  drive  testing.  Figure  35  shows  the  direct  drive  measurement^ 

As  mentioned  above,  this  level  CW 

by  the  filter,  and  does  not  account  for  the 

MOV. 

The  insertion  loss  is  calculated  per  equation 
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dB  (A/Hz) 


FREQUENCY  (Hz) 
EXTRAPOLATED  TO  HEMP 
TP:  B15 


XI  TEST 


FREQUENCY  (Hz) 
EXTRAPOLATED  TO  HEMP 
TP:  B17 


e: 


10< 


Figure  35. 


io» 

FREQUENCY  (Hz) 
A47A 


3.2  DIRECT  DRIVE  OF  THE  TELEPHONE  LINE 


The  insertion  loss  values  calculated  per  equation  (17)  are  listed  in  Table  8. 


3.3  DIRECT  DRIVE  OF  EER  PENETRATIONS 
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SECTION  4 
CONCLUSIONS 


are  less  than 

^■■■■■■■MHHHKeonservative  error  bound  is  12  dB.  These  stresses 
P**J^^5d5ofshield5  and  unshielded  cables,  and  on  wire  bundles.  Note 
that  these  cable  bundles  currents  are  well  below  the  allowed  singlewirestresses  allowed 
in  MTT-STD-188-125.  No  single  wire  currents  were  measured  ^J^mi^Pbecause 

such  meas 


Thus  the  MIL-STD-188-125  field  attmuatignrggui^ents  can  be  met  de- 
the  fact  Th'“  mMns  that 


'4MHI 


t 


I 


.  Since  the  bulk  direct 

“d  ZESTS’  ZZ5ZSBE  the  analysis  and  r«su.ts  are  documented 

separately. 


AH  test  objectives  were  met.  Lessons  learned  and  recommendations  con¬ 
cerning  test  execution  and  test  equipment  improvements  were  provided  m  the  Test 
Director’s  Report  (Reference  l). 


65 


CHAPTER  5 

LIST  OF  REFERENCES 


1.  Thomas  Zwolinski,  “Test  Director’s  Report 
November  1988. 


HEMP  Hardness 


3.  Wolfgang  Bereuter,  “Processing  Direct  Drive  Measurements” , 

to  appear. 


4.  Wolfgang  Bereuter  and  Robert  Racca,  “HEMP  Fields”,  MRC/COS-M-72,  (“un¬ 
published”),  5  December  1988. 


5.  Wolfgang  Bereuter,  “Proposed  Pass/Fail  Criteria  for  MIL- STD- 188-125  CW  Im¬ 
mersion  Test”,  MRC/COS-N-TBD,  February  1989. 


6.  Wolfgang  Bereuter,  “CWMS  Testing  Threat  Relatability” ,  MRC/COS-B-220  31 
October  1986. 


7.  Wolfgang  Bereuter  and  Robert  Racca,  “Processing  of  CWMS  Field  Map  Mea¬ 
surements”,  MRC/COS-N-131,  9  November  1987. 


APPENDIX  A 


HEMP  STRESS  NORMS 

This  appendix  contains  the  numerical  values  for  the  waveform  descriptors  or  norms 
PKI,  PINT,  RIMPI,  and  RAINT  defined  below. 

PKI  =  T |  i(t)  |<  i  ft  |  I(w)  |  du  =  PINT  (Amp)  Peak  and  its  upper  bound 

RIMPI  -  I  *W  I  dt  (Amp  '  #ec)  Rectified  Impulse 

|  di/dt  |<  1  ^jl  |  wl(w)  |  du  =  DIDT  (Amp/sec)  Derivative  Bound 

RAINT  =  {£  Slo  I  !(«)  I2  dM}’  (Amp  V*«)  Root  Action  Integral 

/0  =  10  kHz,  h  =  100  MHz,  T  =  5  #/sec. 
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HEMP  STRESS  NORMS  (continued) 


TP  PHI 

(A) 


PINT  RAINT  RIMPI 

(A)  {A-y/sec)  (A  •  sec) 


A1  6.67E+01  1.18E+02  2.76E-02 

A2  4.86E+01  8.04E+01  1.67E-02 

A3  2.35E+01  6.65E+01  1.55E-02 

A4  2.09E+01  5.93E+01  1.12E-02 

A5  1.30E+01  4.64E+01  8.92E-03 

A6  2.15E+01  5.92E+01  1.50E-02 

A7  1.74E+01  4.55E+01  1.03E-02 

A8  1.95E+01  5.41E+01  9.90E-03 

A10  2.04E+01  6.38E+01  9.95E-03 


4.26E-05 

2.76E-05 

2.73E-05 

2.23E-05 

2.01E-05 

2.58E-05 

1.94E-05 

1.74E-05 

1.86E-05 


All  4.02E+01 
A12  3.15E+01 

A13  2.90E+01 

A14  1.43E+01 

A15  1.61E+01 

A16  1.71E+01 

A17  8.65E+00 

A18  1.06E+01 

A20  1.57E+01 


8.57E+01  9.71E-03 

6.56E+01  7.41E-03 

6.85E+01  9.45E-03 

4.33E+01  6.61E-03 
4.90E+01  7.34E-03 

5.01E+01  6.48E-03 

3.26E+01  5.06E-03 

4.18E+01  7.45E-03 

4.09E+01  5.69E-03 


1.18E-05 

9.83E-06 

1.19E-05 

1.37E-05 

1.31E-05 

1.00E-05 

9.48E-06 

1.31E-05 

1.04E-05 


A21  3.24E+01  6.12E+01 

A22  4.12E+01  7.97E+01 

A23  3.16E+01  6.53E+01 

A24  1.59E+01  5.29E+01 

A25  1.40E+01  4.65E+01 

A26  1.57E+01  5.54E+01 

A30  1.32E+01  3.23E+01 


5.80E-03  6.55E-06 

7.54E-03  7.10E-06 

6.84E-03  6.38E-06 

6.05E-03  9.40E-06 

6.07E-03  8.46E-06 

6.51E-03  6.70E-06 

3.05E-03  3.15E-06 


DIDT 

(A/sec) 


1.02E+10 

7.37E+09 

7.24E+09 

7.29E+09 

5.77E+09 

5.06E+09 

4.55E+09 

6.54E+09 

9.23E+09 

1.52E+10 

1.21E+10 

1.16E+10 

6.87E+09 

8.17E+09 

8.61E+09 

4.87E+09 

6.56E+09 

8.13E+09 

1.40E+10 

1.85E+10 

1.88E+10 

9.93E+09 

8.81E+09 

1.25E+10 

7.83E+09 
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HEMP  STRESS  NORMS  (continued) 


TP 

PKI 

(A) 

PINT 

(A) 

RAINT 

(A’y/aee) 

RIMPI 
(A  •  sec) 

DIDT 

(A/sec) 

Bll 

3.03E-02 

7.21E-02 

2.23E-05 

4.62E-08 

4.92E+06 

B12 

1.61E-02 

5.47E-02 

1.94E-05 

3.38E-08 

4.52E+06 

B13 

2.51E-02 

6.10E-02 

1.52E-05 

2.97E-08 

5.30E+06 

B14 

9.37E-03* 

2.27E-02 

9.58E-06 

2.04E-08* 

3.27E+06 

B15 

3.65E-03 

1.24E-02 

2.82E-06 

5.69E-09 

1.56E+06 

B17 

9.73E-03 

2.15E-02 

1.01E-05 

1.77E-08 

1.36E+06 

B64 

1.02E-02* 

1.90E-02 

3.02E-05 

2.80E-08* 

1.06E+06 

B66 

1.61E-02 

5.82E-02 

2.02E-05 

3.35E-08 

4.52E+06 

B67 

2.46E-02 

6.22E-02 

1.45E-05 

2.91E-08 

5.85E+06 

B68 

3.86E-02 

8.05E-02 

2.35E-05 

4.72E-08 

6.25E+06 

B69 

1.07E-03* 

4.09E-03 

1.48E-06 

1.61E-09* 

5.64E+05 

B70 

1.84E-03 

4.96E-03 

8.53E-07 

1.37E-09 

9.76E+05 

B71 

7.51E-04 

2.40E-03 

5.11E-07 

1.17E-09 

3.21E+05 

B85 

1.54E-02 

4.11E-02 

1.48E-05 

3.54E-08 

2.54E+06 

^30 

4.42E-03 

1.27E-02 

2.75E-06 

8.50E-09 

2.09E+06 

B31 

6.58E-03 

1.55E-02 

3.69E-06 

1.08E-08 

2.62E+06 

B32 

6.49E-03* 

2.11E-02 

6.72E-06 

1.10E-08* 

3.54E+06 

B33 

4.75E-03* 

1.34E-02 

5.15E-06 

3.11E-09* 

2.16E+06 

B34 

3.61E-03 

1.16E-02 

1.99E-06 

2.45E-09 

2.60E+06 

B36 

1.90E-02 

3.34E-02 

4.51E-06 

5.39E-09 

1.04E+07 

B65 

3.14E-02* 

3.07E-02 

4.70E-05 

1.44E-07* 

3.67E+06 

B72 

4.64E-03 

1.56E-02 

3.10E-06 

5.47E-09 

2.68E+06 

B73 

5.01E-03 

2.01E-02 

2.75E-06 

6.37E-09 

4.22E+06 

B74 

4.28E-03 

1.67E-02 

3.14E-06 

8.44E-09 

3.44E+06 

B75 

1.85E-03 

3.54E-03 

1.18E-06 

4.14E-09 

9.20E+05 

B76 

1.84E-04 

4.86E-04 

6.75E-08 

1.19E-10 

1.22E+05 

B77 

2.91E-03 

6.43E-03 

1.27E-06 

1.56E-09 

2.03E+06 

B86 

2.55E-03* 

7.99E-03 

8.56E-06 

6.86E-09* 

8.01E+05 
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HEMP  STRESS  NORMS  (concluded) 


TP 


PKI 

(A) 


PINT 

(A) 


B49 

2.76E-03 

8.41E-03 

B50 

5.72E-03 

1.40E-02 

B51 

8.46E-03* 

1.48E-02 

B55 

3.58E-03 

1.48E-02 

B62 

5.87E-02 

1.27E-01 

B63 

2.46E-02 

5.74E-02 

B78 

7.73E-03 

1.61E-02 

B79 

6.41E-03* 

1.85E-02 

B80 

4.55E-03 

1.32E-02 

B81 

2.45E-03* 

6.46E-03 

B82 

7.36E-04* 

1.67E-03 

B83 

5.08E-03 

8.38E-03 

B84 

3.12E-03* 

6.21E-03 

G3 

1.91E-02 

4.98E-02 

G8 

3.42E-02* 

2.25E-02 

G9 

4.51E-03* 

1.13E-02 

G10 

8.61E-03* 

1.06E-02 

Gil 

1.15E-03* 

2.98E-03 

G12 

6.40E-03* 

1.21E-02 

G13 

4.07E-03 

1.22E-02 

G14 

9.54E-02* 

4.83E-02 

G15 

7.63E-03* 

1.32E-02 

G16 

6.67E-03* 

1.60E-02 

G17 

1.03E-02* 

2.37E-02 

G18 

1.18E-02 

2.58E-02 

G19 

1.03E-02* 

1.95E-02 

G20 

5.87E-03* 

1.34E-02 

G21 

4.17E-03* 

9.98E-03 

G31 

1.68E-02 

4.13E-02 

G32 

1.65E-02* 

2.97E-02 

G33 

1.32E-02* 

1.84E-02 

G34 

1.10E-02* 

1.24E-02 

"  W100 

2.30E-03* 

2.87E-03 

W101 

1.86E-03* 

1.79E-03 

raint  rimpi 

(A-y/sec)  (A  •  sec) 


1.19E-06 

1.92E-09 

3.40E-06 

8.55E-09 

6.05E-06 

1.76E-08* 

2.31E-06 

5.38E-09 

1.47E-05 

1.34E-08 

7.75E-06 

1.01E-08 

2.95E-06 

4.72E-09 

6.07E-06 

3.63E-09* 

1.62E-06 

2.63E-09 

4.63E-06 

6.67E-09* 

5.43E-07 

1.18E-09* 

3.32E-06 

1.12E-08 

3.45E-06 

1.06E-08* 

7.59E-06 

9.25E-09 

3.95E-05 

1.41E-07* 

4.77E-06 

•  1.35E-08* 

1.05E-05 

3.60E-08* 

1.35E-06 

1.51E-09* 

6.55E-06 

1.94E-08* 

2. 64  E- 06 

6.35E-09 

1.36E-04 

4.62E-07* 

1.63E-05 

2.97E-08* 

4.91E-06 

1.59E-08* 

1.58E-05 

3.45E-08* 

6.93E-06 

2.07E-08 

9.01E-06 

2.65E-08* 

5.17E-06 

1.50E-08* 

3.88E-06 

1.16E-08* 

8.82E-06 

1.10E-08 

1.24E-05 

4.1  IE-08* 

1.68E-05 

4.45E-08* 

1.32E-05 

4.50E-08* 

5.89E-06 

1.05E-08* 

1.63E-06 

5.26E-09* 

POOR  FOURIER  INVERSE  TRANSFORM  (FTE  >  .3) 


DIDT 

(A/sec) 


1.62E+06 

2.52E+06 

2.36E+06 

2.60E+06 

3.28E+07 

1.80E+07 

3.64E+06 

3.38E+06 

2.70E+06 

1.50E+06 

3.37E+05 

2.03E+06 

1.08E+06 

2.22E+07 

1.84E+06 

2.88E+06 

2.30E+06 

7.99E+05 

1.97E+06 

3.14E+06 

2.47E+06 

2.36E+06 

4.25E+06 

3.90E+06 

6.40E+06 

3.49E+06 

2.29E+06 

1.60E+06 

1.41E+07 

8.77E+06 

1.69E+06 

2.53E+06 

1.04E+05 

1.08E+05 
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